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' We propose a new criterium for saturation of the density of partons both in nucleons and 

' , nuclei. It is applicable to any multiple scattering model which would be used to compute the 

Qh. 

^ ■ number of strings exchanged in ep and eA collisions. The criterium is based on percolation 

of strings, and the onset of percolation is estimated from expectations coming from the 
• 1—1 , 

^ ' study of heavy ion collisions at high energies. We interpret this onset as an indication of 

H : 

saturation of the density of partons in the wave function of the hadron. In order to produce 
quantitative results, a particular model fitted to describe present HERA data and generalized 
to the nuclear case is used. Nevertheless, with the number of scatterings controlled by the 
relation between inclusive and diffractive processes, conclusions are weakly model-dependent 
as long as different models are tuned to describe the experimental data. This constitutes a 
new approach, based on the eikonal description of soft hadronic collisions, and different from 
others which employ either perturbative QCD ideas or semiclassical methods. It offers an 
alternative picture for saturation in the small region. 
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Much interest has recently been devoted to the saturation of partonic densities 
i.e. the change in the increase of the partonic densities from power-hke to logarithmic 
or constant with decreasing parton momentum fraction x, both in nucleons [Q, ^ 
and in nuclei §|, 01- From the point of view of experimental data on lepton-hadron 
scattering, the most striking feature was the change in the logarithmic slope of the 
proton structure function dF2/d\n {Q"^) at x ~ 10~^, the so-called Caldwell plot (now 
known to be mainly due to a Q"^ — x correlation), but the situation is not conclusive: 
Nucleon data can be described not only in approaches which consider saturation 0, ^, 

but also satisfactorily accommodated in the usual global fits (also available for 
nuclei |T^) that consider the standard QCD evolution or resummation [^, starting 
from initial conditions at low photon virtualities which do not include saturation 
(see |T2| for an application to the Caldwell plot and [^] for a discussion on the present 
situation) . 

From a theoretical point of view, the saturation regime is a very interesting one 
characterized by a small coupling constant and high occupation numbers, where a 
semiclassical description in terms of fields has been proposed 0. Different models of- 
fer explanations based on multiple scattering (i.e. unitarization) or gluon interaction, 
both in the case of nucleons H, 0, ^ and nuclei [jl], |^, |], 0. These two approxima- 
tions to the problem are equivalent (see e.g. [|l^) in different reference frames, but the 
models predict the onset of saturation in different kinematical regions and the satura- 
tion features are also diverse. In this short note we will essay another approach to the 
problem, inherited from multiparticle production in nucleus-nucleus {AB) collisions at 
high energies and applicable to any model formulated in terms of multiple scatterings. 

The concept of saturation, not of the density of partons in the hadronic wave func- 
tion but of the number of partons produced in the collision, was proposed some time 



ago [|I^ in AB collisions at high energies and has been reconsidered recently in the 
context of the search of the Quark Gluon Plasma (QGP); such high partonic density 
should provide the initial condition for the possible thermalization of the created sys- 
tem. Several related ideas have been used to compute the multiplicity of produced 
particles in AB collisions at the Relativistic Heavy Ion Collider (RHIC) at BNL and 
at the future Large Hadron Collider (LHC) at CERN [p!6| , For example, in p!8[ 
perturbative QCD (pQCD) is used to compute the initial number of gluons, quarks 
and antiquarks, which are limited according to the simple geometrical criterium that 
the number of partons per unit of transverse space times their transverse dimension 
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(oc cannot be greater than 1 (see for other attempts in this direction). Be- 

sides, the semiclassical methods used in have also been employed to estimate the 
initial number of gluons in a heavy ion collision [ 2D | . 

On the other hand and in the framework of string models for soft multiparticle pro- 
duction (see and references therein), a simple geometrical criterium for saturation 
has been proposed. In these models particle production comes from string breaking, 
strings which are considered, in a first approximation, as formed and decayed inde- 
pendently. As the number of binary nucleon-nucleon collisions (each one producing 
2 strings [^, |22[) increases with increasing centrality, energy or nuclear mass, this 
approximation should break down. The onset of this phenomenon can be estimated 
considering strings with a certain area in the transverse space of the collision, and 
taking into account the possibility of two-dimensional percolation of the strings when 
they overlap in this transverse space. Percolation is a second order phase transition 
which takes place when clusters of overlapping strings, with a size of the order of the 
total transverse area available, appear. This idea has been proposed in AB collisions 
23[| and applied to signatures of QGP . 

The purpose of this note is to use percolation of strings as an indication for the 
onset of saturation of the density of partons in nucleons and nuclei, a quantity which 
in our case is not directly related with the partonic densities measured in DIS, as such 
identification ||, Q can only be done at high Q"^, and our approach will be devoted 
to the low regime. For this, we need a multiple exchange model for ep collisions 
which allows us to compute the number of binary collisions, generalize it to the nuclear 
case, translate the number of collisions to a number of strings and estimate the density 
of strings to compute whether percolation takes place or not. The method can be 
applied to any multiple scattering model, and the results in any of these models should 
be quite the same (within the uncertainties due to the extrapolation of the model to 
nuclei and to higher energies or smaller x) as long as the model is able to describe the 
fully inclusive and diffractive experimental data on ep collisions, see comments below. 

Let us give a brief description of the model developed in 0] , which is the particular 
one we are going to use to compute the number of binary collisions and then of strings 
in order to give quantitative predictions. The goal of the model was the description of 
total and diffractive data on ep scattering at low and moderate and small x. This 
region is where unitarity corrections are more important and where the transition from 
non-perturbative to perturbative QCD takes place. 
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In the proton rest frame, the virtual photon coming from the lepton fluctuates into 
a qq state. Then this hadronic state interacts with the proton. Unitarity corrections 
are described by the multiple scattering of the qq fluctuation with the proton. In a 
quasieikonal approach the total cross section is given by 

^L'(s, Q^) = MQ^)J ^i^^Mzg(£i«LMJ , (1) 

where giQ"^) is the 7* — {qq) coupling, 2x(s,Q^,6) is the elementary [qq) —p cross 
section at fixed impact parameter b, and C = 1.5 is a parameter taking into account 
the diffractive dissociation of the proton. For small sizes r of the qq pair, x oc from 
pQCD calculations. As oc for these small sizes x 0^ VQ^- For large sizes of 

the fluctuation no Q^-dependence is expected. In [|| two components, corresponding 
to small (S) and large (L) sizes of the qq pair, were taken into account, x = Xl + Xs- 
The fact that Xs while xl does not depend on makes the correction terms 

in Eq. (Q) more important for the L part than for the 5* one. So, more scatterings are 
present - in average - for the L than for the S component; for this reason and also 
due to the fact that we will consider small Q^, only the L component, which is the 
dominant one ^ for < 2 (GeV/c)^, will be used in the actual computations. The 
energy dependence of these x's is given by a single pomeron of intercept A = 0.2, 

Here, ^ = In^^, \l = Rl + a'pi with i?i = 3 GeV-^, and a'p = 0.25 GeV-^, the 
slope of the pomeron trajectory, gives the Ins behavior of the total cross section for 
very large s; besides, Cl = 0.56 GeV^^ and Sq = 0.79 GeV^. The variable ^ is chosen 
so that X for large ^ sq and x oc (s/sq)^ for 0; in this way, the model 

can be used for photoproduction. For the 5* part, similar expressions were used in |^ 
with an extra factor in Eq. (^. 

The description of diffraction is a very important ingredient of the model. It is given 
by quadratic and higher order terms in x in the expansion of Eq. ([^). Thus, the ratio 
adiff/dtot controls the unitarity (multiple scattering) corrections, i.e. the number of 
scatterings and strings; this idea has been used to compute nuclear structure functions 
from a description of diffraction at HERA, see e.g. |2^. So, any multiple scattering 
model able to reproduce the experimental data on this ratio should produce roughly 
the same number of scatterings (strings) and, consequently, the same predictions for 
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the onset of percolation and saturation. A triple pomeron term was introduced in 0] 
in order to reproduce large mass diffractive processes. This term is another source of 
shadowing corrections to the total cross section. It will be used in the actual computa- 
tions, see 0] for the full expressions and parameters. The reggeon contribution which 
appears in decreases with increasing energy and is negligible at the energies under 
consideration, so it has been ignored. 

The model in |Q has 9 free parameters that were fitted to experimental data on 
diffractive and total ep cross sections for < < 10 (GeV/c)^ and 10"^ ^ 2; ^ 10~^. 
Once the parameters of the model are fitted, it is possible to know the mean number 



of collisions p2| : 



n = } ^ n 



Jd% 2Cx{s,Q^b) 
Jd%'[l-exp{-2Cx{s,Q^b')}] 



(3) 



where, for n > 1, 



(T„(s, Q M = —F, j exp(-2(:7x). (4) 

Notice that in these expressions, ^(s, Q^, &) contains the triple pomeron contribution, 
so cuts in different branchings of one single fan diagram (i.e. one single tree of triple 
pomeron couplings) are included in the same (T„, which thus corresponds to the ex- 
change of n fan diagrams, each of them cut in one o more than one of its branches. 
Thus Eq. (1) is a conservative estimation, as these cuts could give rise to a larger 
number of strings. Besides, all our expressions are asymptotic ones, not considering 
energy-momentum conservation (which could reduce slightly the number of collisions 
at the lowest energies). 

Neglecting isospin at high energies, the generalization of any multiple scattering 
model formulated for ep to the case of eA collisions is straightforward in the Glauber- 
Gribov approach |^6|]: The number of qq — nucleon collisions (the number of partici- 
pating nucleons of A) in this case, is given [27] in terms of the inelastic non-diffractive 
cross sections by 

M = A af/at^ (X (5) 

with 

= 9{Q') J d'b (1 - exp {-AT^{b)aj:'/g{Q')}) , (6) 
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TA{b) = J^^dzpA{z,b) the profile function normalized to 1 taken from ||2^ and 

CO „ _ 



n=l • 

gm 

c 



Jd%[l-exp{-2Cxis,Q^b)}]. (7) 



So, the total number of collisions is given by {ricoii) = {npart) n. As previously com- 
mented, in the actual computations we will only use the L component, Eq. (^, as it 
is the dominant one |Q for < 2 (GeV/c)^ where our calculations will be done. 

At this point, it could be argued that using the model in (or any other multiple 
scattering model) there is a possibility to study saturation of the density of partons, 
both in nucleons and in nuclei, simply looking at the point in which amplitudes in im- 
pact parameter space become energy independent, or alternatively the point in which 
cross sections reach a regime in which their energy behavior becomes identical to that 
of the size of the target (expanding in the case of a nucleon). Nevertheless, the gener- 
alization of 1^ to nuclei is not so obvious: ours is a very simple one, but more rigorous 
generalizations ||29| also rely in simplifications of the exact Gribov calculus pO] or 



Glauber-Gribov theory |26|. So we think that an estimation as the one we perform, 
based on geometrical criteria, is worthy, of simple and general applicability, and may 
provide, as in the case of nucleus-nucleus collisions, an indication of the onset of a high 
density, non-linear regime. 

Let us establish now our criterium for saturation of the density of partons in the 
wave function of the target. As it was said, percolation is a non-thermal second order 
phase transition, which takes place when clusters of overlapping objects acquire a size 



comparable to the total size available [31|. In our case the space is the transverse 



dimension available for the collision, and the overlapping objects are strings. The 
parameter which controls the onset of percolation is the dimensionless string density 

r] = Nt/T (8) 

(which may be related [O with the dimensionless density of gluons found in semiclas- 



sical models), with 

= 2 {n,oii) (9) 

the number of strings exchanged in the collision (each collision gives rise to two strings 
due to the pomeron dominance at high energies, see [^, |22[), t = hTq is the transverse 
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dimension of the string^, with tq — 0.20 -i- 0.25 fm as extracted from phenomenology 



23| , |2^ , and T the total transverse area available for the collision. This last quantity 
is not known and could depend on the virtuality of the fluctuation Q^; however, for 
small and moderate Q^, it can be estimated to be the typical size of the vector meson 
in which the virtual photon fluctuates (as this is the smaller object in the interaction). 



So, we will use T = 1 fm^ (a radius yT/n ~ 0.56 fm). Also, a size varying with the 
energy, in the spirit of an expanding proton, could be explored; for example, a size 
increasing with increasing energy would slow the corresponding increase of the density 
of strings but, for simplification, we will employ a fixed size. 

The critical value for rj where percolation takes place, has been computed using 
different methods and depends quite strongly on the profile of the nucleus (i.e. on 
the distribution of the overlapping objects inside the available transverse space). For 



continuum two-dimensional percolation and from we take 77^ ~ 1.12 1.50. 

Defining the string density as n = N/ T and allowing for the different values of ro and 
r/c, we find a critical string density 

ric — Q -T- 12 strings/fm^. (10) 

With this critical string density, it is tempting to estimate the behavior of the at 
which, for a fixed s, saturation of the density of partons takes place in this approach, 
Q'sat [S 05 S H) B 0) mi- However, in our case, being xl almost independent on Q^, 
the density of strings is also almost Q^-independent for fixed s; besides, the model is 
only valid for small and has not been designed for Q^-evolution. More significant 
in our model is the value of x or s where, for fixed and A, saturation takes place, 
Xsat or Ssat respectively. Considering neither triple pomeron nor reggeon contributions 
and approximating in Eq. (|^) a^^"^ ^ g{Q^) t^RqA^/^, with _Ro — 1-2 fm, it is found 
that x,at oc [QV(^o + Q')] A^l^^^') and s^at oc (sq + Q^) A-i/{3A)^ i/(3A) = 5/3. So, for 

^ 0, Xsat increases linearly with increasing while Sgat is roughly constant; these 
qualitative features will be observed in the numerical results. 

Let us turn to numerical evaluations. Using Eqs. (^-(0), d^) and (0), we can now 
compute the string density for different hadronic targets, Q^, and x or s. When this 
density becomes larger than the critical value, Eq. (p!oD , percolation will take place, 
which we will interpret as a signal of the onset of saturation of the density of partons 

^In this approach this transverse dimension plays the role of an intrinsic scale of soft physics, 

— > 0, where a description in terms of pQCD degrees of freedom becomes dubious. 
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in the target. In Fig. 1 we present results for the string density in 7* — A colhsions 
versus x for = 0.1 and 1 (GeV/c)^, with A = 1, 9, 56 and 207 (corresponding to 
p, Be, Fe and Pb respectively). In Fig. 2 the same quantity is presented versus s for 

= and 2 (GeV/c)^. Some comments are in order: First, from Fig. 1 it looks 
as if saturation (percolation) is favored by a higher Q^, apparently in contradiction 
to what is commonly expected. This is due to the fact that, as we have said, the 
number of rescatterings is hardly dependent on - as can also be observed in Fig. 
2, and that s is the variable which controls this number (indeed, in Fig. 2 it can be 
seen that Sgat increases slightly with increasing Q^, as expected). So, for a fixed s at 
which percolation occurs, the higher the the higher the Xgat (in agreement with the 
naive expectations in the previous paragraph). Second, the dependence of Ssat on A, 
parametrized as A~°', is found to be stronger in the numerical computations (a ^ 5/2) 
than the power 5/3 estimated in the previous paragraph. This discrepancy is due to 
the triple pomeron contribution included in the numerical computations, which appear 
as a denominator in Eq. (H), diminishing the 'effective' A which appears in Ssat and 
thus making a = 1/(3 A) larger. 

To summarize, a criterium for saturation in the small region applicable to any 
multiple scattering model, has been presented. To produce quantitative results, a mul- 
tiple scattering model for 7* — p collisions in this region 0] has been generalized 
to the nuclear case, and used to compute the number of exchanged strings. As mul- 
tiple scattering (and thus the number of produced strings) is controlled by the ratio 
a°'*-^-^/cr*°*, this number is related to experimental data on diffraction and the actual 
realization of the model is not crucial to compute the string densities as long as it 
reproduces the experimental data. Employing the ideas of percolation of strings taken 
from Heavy Ion Physics, the kinematical regions for the onset of percolation, which 
has been interpreted as saturation of the density of partons in the target^, have been 
calculated. This constitutes a new approach, based on Regge phenomenology and dif- 
ferent from others which use either pQCD ideas or semiclassical methods; it offers an 
alternative picture, based in hadronic degrees of freedom (strings), for saturation in 
the small region. In view of the results presented in the Figures for the onset of 
percolation, which for large nuclei may appear at not so small x, saturation could be 



observed in future eA colliders |33], and the effects of this second order phase transition 



^As commented previously, we are working in the low regime, so this quantity cannot be 
identified with the partonic densities measured in DIS. 
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visible in correlations (as proposed in AB collisions ||34|| ). 
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Figure captions: 



Fig. 1. String density in 7* — p, Be, Fe and Pb collisions versus x for = 0.1 
(GeV/c)^ (solid lines) and = 1 (GeV/c)^ (dashed lines). Dotted lines are the 
bounds on the critical string density for percolation of strings. 

Fig. 2. String density in 7* — Be, Fe and Pb collisions versus s for photoproduction 
(solid lines) and — 2 (GeV/c)^ (dashed lines). Dotted lines are the bounds on the 
critical string density for percolation of strings. 
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